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ABSTRACT 


A formula for transformation of solid spherical harmonics under 
coordinate changes makes it possible to convert potential functions, 
expressed m a series of these harmonics, to representations valid 
in different coordinate systems 

Such a formula has several geophysical applications when the po- 
tential of a celestial body has been determined relative to one parti- 
cular set of coordinates A coordinate change may improve the con- 
vergence of the series m a given region, or such a change may be con- 
venient for other reasons 

There are also applications to micro -gravitational interactions 
among non-spherical, but rotationally symmetric bodies, allowing the 
analytic calculation of quantities which depend on potential functions, 
such as gravitational forces and torques 

In this work the necessary transformation formulas are developed 
and applied to the task of analytically calculating the force between 
two homogeneous hemispheres Two special cases are worked out m de- 
tail as examples where the coordinate changes involve both translation 
and rotation 


Thesis Supervisor. Stephen J Madden, Jr 
Title Lecturer 
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CHAPTER 1 
INTRODUCTION 


1 1 Motivation 

The motivation for this research was the need for a satisfactory 
method of calculating forces and torques caused by gravitational at- 
traction among bodies on the laboratory scale 

Historically, the calculation of gravitational forces has been 
limited to relatively large forces - those involving celestial bodies, 
where all, or all but one of the masses, can effectively be treated 
as a point mass Under these conditions , the calculations reduce to 
a simple gradient operation on the scalar potential. 

Laboratory experiments which measure the magnitude of gravitational 
forces and torques using test masses, such as the Cavendish experiment, 
have in the past involved at least one spherical test body of uniform 
density The problem again reduces to one of determining the potential 
and its gradient 

Currently, however, attempts are being made to greatly improve 
the accuracy of the Cavendish measurement m order to make a rela- 
tivistic check of the dependence of the gravitational constant on the 
potential The use of non-sphencal test masses appears to be one 
way m which the sensitivity can be improved (Lee, 1970) Calculating 
the forces and torques in this case involves a triple integral, or the 
case of homogenous mass , a double integral over the volume or surface , 
respectively, of the passive mass. 
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V V dx 


F 


/ 

volume 


=* J V n dS, (1-1) 

surface 

where is the mass density, n is the outward normal to the surface, 
and V is the potential due to the gravitating mass(es) In general 
this calculation is formidable - difficult to carry out analytically, 
and computer approximation schemes do not efficiently yield the 
accuracy required for a relativistic measurement* 

1 2 Objectives 

The purpose of this investigation is to develop an analytical 
method for calculating the gravitational interactions among rotationally 
symmetric bodies of uniform density, permitting the evaluation of the 
forces to the desired accuracy Other quantities which depend on the 
potential, such as torques and partial derivatives thereof, can be 
treated m the same manner 

1 3 Method 

To facilitate the evaluation of the surface integral m equation 
(1-1) , the coordinate system is oriented to take full advantage of the 
symmetry of the passive body The coordinates and coordinate system 
originally selected are best suited for expression of the potential 
function of the gravitating mass(es). Once the potential has been 
determined, the reference frame will, m general, have to be changed 
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before permitting the integrations. Although the scalar potential re- 
mains unchanged under such a transformation, individual terms m a 
series expansion of the potential will change with the coordinate trans- 
formation 

In this thesis a potential expansion xn solid spherical harmonics 
is transformed under translation and rotation of coordinates, permitting 
the relatively easy evaluation of the surface integrals The trans- 
formation will involve one translation and one rotation, unless the 
translation is limited to the direction of the axis of symmetry 
of the gravitating mass, m which case two rotations are required in 
addition to the special translation* 

The transformation of the expansion requires the transformation 
of the coefficients in the series, which m turn necessitates the 
transformation of the solid spherical harmonics under the coordinate 
change 

1 4 Other Applications 

The transformation of a series expansion in solid spherical har- 
monics appears to have several applications other than the one des- 
cribed in 1.1 The potentials of the earth and moon, for example, 
which are originally determined m terms of a particular coordinate 
system, can then be written m any other coordinate system desired 
This may facilitate its evaluation m some regions, either because 
the convergence is not satisfactory, or to capitalize on a particular 
geometrical configuration 


11 



PRECEDING PAGE BLANK NOT FILM1 


CHAPTER 2 

TRANSFORMATION OF SOLID SPHERICAL 
HARMONICS DUE TO TRANSLATION OF AXES 

2 1 Coordinate Systems 

A pure translation of the origin of the x, y, z coordinate 
system (Figure 2-1) to a new system (£,u,£) 1S represented by a trans- 
lation vector £ A field point P is located by the position vector 
r m the old system and a vector R m the new system 



Figure 2-1 Relationship Between Old and New Coordinate Systems 

In spherical coordinates, the angles 8, X and ct are the co- 
latitudes, (j>, y and 3 the longitudes, and r, p and R the magnitudes 
of the three vectors r, £ and R, respectively The coordinates of 
r are indicated m Figure 2-2 as an example 
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Figure 2-2 Definition of Spherical Coordinates 


Thus , 


r “ R + £ , 


where 


r = 



r sm 0 
r sm 0 
r cos 


" 5 " 


R sin a 

T1 

= 

R sm a 

K 


R cos 

- _ 


— 


P 


r 

p sm X cos y 
p sm X sin y 
p cos X 


cos 4> 
sm <{> 

e 

cos 8 
sm 8 
a 




and 


0 < r, R, p 


0 < 0, a, A < it 


0 < $ , 3 , y ^ 2 tt 


2.2 Solid Spherical Harmonics 

A solid spherical harmonic is defined as any homogeneous func- 
tion f(x,y,z) of degree n which also satisfies Laplace’s equation, 

V 2 f = 0 (2-2) 

If harmonics of the form 

r^P^(c°s 6)e im< ^ n=0,l,2, 

m=-n , -n+1 , ,n-l,n 

are considered, then p must assume one of two values, 
p = n, -n-1 , 

since P™(cos 0)e lin<i> satisfies equation (2-2) (MacRobert, 1947, p 150) 
The P^'s are known as associated Legendre polynomials, defined by 

pj(v) = (-l) m (l-v 2 ) m/2 P n (v) , (2-3) 
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where the P n *s are the familiar Legendre polynomials They can be 
determined from Rodrigues 1 formula 


V”> 


2V 


dv 


n 


(v 2 -l) n . 


Several important relationships among the P^'s are noted here for 
later use The first relates a polynomial of negative order to the 
corresponding polynomial of positive order (Jackson, 1966), 



(v) 


(-l) m (n-m) ' 
(n+m) 1 


P m (v) 

n 


m > 0 


(2-4) 


Note that m the special case m = 0, 


p”(vl - p n (») 

Other useful expressions are (Abramowitz , 1965) 

pjj(v) = (-I) n+ V(-V) (2-5) 

p m = (-l) lll (2iii) 1 (1 _ y 2 m/2 

m m'2 m 

P m+1 < v) = < 2m+1 > v PjN) , (2-6) 

and the recursive relations 

(n-m) P™ (v) - (2n-l) v p ™_! (v) + (n+m-l)P™_ 2 (v) = 0 (2-7) 

(n-m) (n+m+l)P™(v)+2(m+l)v(l-v 2 ) _1 ' /2 p™ +1 (v)+p™ +2 (v) = 0 

(2-8) 
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Several P n (v) f s and P^{v)*s are evaluated for later reference. 


P 0 (v> 

= 1 

P^(v) = -(l-v 2 ) 1/2 

P X (V) 

= V 

P 3 (V) = -3(l-v 2 ) 1/2 v 


- |(3V 2 -1) 

P 2 (v) = 3 (1-v 2 ) 

P 3 (V) 

= J{5v 3 -3v) 

P 3 (v) = - |(l-v 2 ) 1/2 

P 4 (V) 

- i(35v 4 -30v 2 +3> 

P 2 (v) = 15 (1-v 2 ) v 

P 5 (v) 

= -|(63v 5 -70v 3 +15v) 

p|(v) = -15 (1-v 2 ) 3 / 2 


A well known and useful generating function for harmonics involving 
the P n ! s (Erdelyi , v Q i 1, 1953, n J.S4) is 


(l+h 2 -2vh ) 1/2 


n=0 


l h- n_1 p 

n=0 


n 


(v) 


| h | < smaller of 

|v+(v 2 -l) 1/2 | 

|h| > larger of 

| v+{v 2 -l) 1/2 | / 
(2-9) 
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and finally, the addition theorem of spherical harmonics, involving 
the angle between two vectors r and r 1 (Jackson, 1966, p.67)# I s 

P n (cos 1 J 1 ) - P n (cose) (cos 9 1 ) 


+2 


nv 


Y fo-rc) 1 
(n+rtiy ' 


P®(cos 


0)p“(COS 


0 ’ ) cos m($-$ f ) 


(2-10) 


where 


cos i|> = cos 0 cos 0' + sm 0 sm 0* cos { 4>— 4 1 1 ) 

2 3 Transformation of Harmonics Involving Powers of r n 

In order to transform the coefficients in solid spherical harmonic 
expansions with powers of r 11 and r n ^ , the spherical harmonics them- 
selves must be transformed m the opposite sense for each case That 
is, the old harmonics must be found m terms of the new ones 

If the method developed by Aardoom (Aardoom, 1969) is adopted, use is 
made of the finite binomial expansion 

(u T r) n = (u T R+u T p) n = l (P) (u T R) n ' k (u T p) k , (2-11) 

k=0 K 


where 


= n' 

V " (n-k) ’k’ ' 
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— T 

and (u r) is the Hermitian scalar product of a special dummy complex 
vector u, and r* The bar over u indicates its complex conjugate and 
the superscript T stands for its transpose 

The vector u is defined m terms of a free parameter t, 

t ^ 0 , 


u = 


1-t* 

i(i+t 2 > 

-2t 


and the Hermitian scalar products can be identified with the generating 
function in equation (A-ll) of Appendix A Thus, 


(u T r) n = t n | ll) H {r , cos 6,$)t m , 

m=-n n,m 


(1) 


(2-12) 


where 


(1), T /_ a ^ . (-l) In (-2) n n , „n„m, 

H n,m (r ' cos = (n+mp* r P n (cos 6)e * (2 13) 


Using equation (2-12) in three places m (2-11) , and dividing by 


t n . 


i {1) \ = i (j) "r ll, H n . k hiRit 11 1 ll \ (p)^ 

m=-n n ' m k=0 K h=-n+k n K ' n 3 =~k 


n _ n-k 


(i), 


uh v (1), 


(2-14) 


With a redefinition of indices, (2-14) can be written as 


5 r, 5 r, „ l -k+n+p 


l H n ‘ E t p [ (£) l H k H n-k D -£ (p) ' 

i=-n n ' m p=-n k=0 k . _ -k k ' 4 n K ' p 1 

^ »lk+p-n 


(2-15) 
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where 


si and l|* 

<b * b 

are taken as the smaller and the larger, respectively, of the two 
quantities (a,b) beside the bracket. 

Since t is a free parameter, the coefficients of t m in equation 
(2-15) can be equated In spherical coordinates, making use of equa- 
tion (2-13) , the result is 


r np ™(cos 9)e 


imc|> 


n (-k+n+m,.. . „ , . 

I I , (1) K^ r i(p,cos X, Y )R k p£(cos a)e lAp , 
k=0 


n,m ' 


(k+m-n 


(2-16) 


where 


(l) K k,A 
n ,m 


(p / cos X,Y) 


r n+nw n-k_m-£ 
C k+Jt )p P n-k 


(cos X)e l(ra " i)Y 


(2-17) 


k Z 

Note that xs a function only of the transformation variables p, X 

and y 


2 4 Transformation of Harmonics Involving Powers of r~ n 1 

The transformation of the harmonics involving inverse powers of 
r is treated in a manner analogous to that used in section (2 3), 
beginning this time with an infinite binomial expansion. 


(u T r) _n_1 


(u T R+u T p) n 1 


{ I (-l) k ( n J k )(5 T R)- n - k - 1 (G T p) k 

k=0 n 

i (-i) k ( n ; k )(u T p)" n *- k - i ( a T R) k 

k=o n 




( 2 - 18 ) 
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Thus r the series takes one of two different forms, whichever will con- 
verge. The convergence criteria themselves are difficult to evaluate 
due to the involvement of the dummy vector u In the special case of 
a translation along the z axis (X=0 ) , it can be shown that the cri- 
teria reduce to the question of £ being less than or greater than one. 
This has not been directly shown to be true in general, except that 
an independent derivation by Madden (Madden, 1971) of the transforma- 
tion of the harmonics yields the same results as will be obtained 
here, with the criteria 

p > R 

holding in general Hence, these convergence criteria are adopted 
The left-hand side of (2-18) can be identified with (A-13) , so 

that 



-n-1 


-n-1 


n 

I 

m=-n 


( 2 ) 


H (r,cos 0 , 4>) t 

Ti f m 


m 


(2-19) 


where 

(2 V (r,cos 0,*) = r' n 'V(cos eje 1 " 1 '*’ (2-20) 

n,m (-2) n ill n 

ntl 

Making use of (2-12) and (2-19) and multiplying by t , equation 
(2-18) can be expressed as 


n 

l 


m=-n 




- i (-i> k (”; k ) 

k=0 n 


“f (2)„ fE 

"+ k - h 


D =-k 


( 1 ) 


H. 




£ 

i R I 


( 2 - 21 ) 
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where 



Changing indices and equating coefficients of t m on each side of (2-21) 
gives 


l2) »„ - I <-i> k o T (2, H„ +k . 

n ' m “ k=0 n *4-k a+k ^ 

r n 

R 

£ 

d)„ 

P. 

R 


» 


. ^ 


Introducing spherical coordinates and expressions of the form of (2-13) 
and (2-20) yields 


.,-n-k-l i 


R 


l 

|_m-£ 

k+n 

k 


(cos a)e 


• [m-l] 1 


( 2 - 22 ) 


where 



f n k 1 


m-JT 


( 2 ) --k f it „ * j 

n,m ~ 1 K n-m ; 

P 

-n-k-1 

P 

_ 1 
‘ k = 

1 1 . lY 

(cos X)e , {2 _ 23) 


L P J 


k+n_ 



and 



p < R 
p > R 


In the special case when the translation is in the positive z 
direction (1=0), equations (2-16) and (2-22) reduce to ; 
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(2-24) 


r n P™{cos 8 ) e 


im$ _ 



n-k R k p m 


, . uaf5 

(cos a)e 


t 


and 


l {-l) k ( k+ ^" It, )P k R" n " k “ 1 Pn+k (cos “ )e 

k=0 K 


im0 

p < R (2-25) 


r - n — V“( cos 6)e im * = 


l ( - 1 ) k+m ( ) P " n_k " 1 R k P J ( cos «)e 

k=0 k+m 


im8 


p > R (2-26) 
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CHAPTER 3 

TRANSFORMATION OF SPHERICAL 
HARMONICS DUE TO ROTATION OF 
THE COORDINATE AXES 


3 1 The Rotation Matrix 

To effect a transformation of coefficients due to a rotation of 
coordinate axes (Chapter 4) , the new harmonics are first deter- 
mined m terms of the old ones A particular rotation of axes defines 
a transformation matrix, M, such that 

R(R,a,(3) = M r (r , 0 , <{>) , (3-1) 

where M is formed by successive rotations through three Euler 
angles 1) ft about the z axis, 2) I about the newer x axis, and 
3) to about the new z axis, all m a right-handed sense (Figure 3-1) 



0 < ft £ 7T 

0 i 1 1 tf/2 

0 4 < to < 2 it 


Figure 3-1 Euler Angles Associated with a Rotation of Coordinate Axes 
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Adopting the method of Courant and Hilbert {Courant and Hilbert, 
Vol 1, pp 535-545), a skew symmetric matrix. A, is constructed in 
terms of three parameters, q^, q 2 and q^. 


- r 


0 


q 3 


q 2 


A = 


" q 3 


0 


q l 


L q 2 ~ q l 


(3-2) 


Three additional parameters, q 4 , w and v, are introduced by 


, 2 
w = (q 1 

* 4 ♦ 4>v* 

(3-3) 

, 2 
v = ( q i 

*4*4* *4> 1/2 - 

(3-4) 


where q^, q 2 » q 3 and q 4 are all real. The immediate objective is to 
show that the orthogonal matrix M can be expressed as 

M = (q 4 I M + A) (q 4 I M _- A)" 1 , (3-5) 

* 

the exponent of minus one indicating the inverse of the matrix, and Ij 
representing the identity matrix 

Note that (q 4 I M +A) and commute Hence, (q^I^+A) and 

(q 4 Ij 4 “A) commute, since, in general, D 1 C = CD 1 , if CD = DC 
It follows that 

M T = (q 4 I M - A T )" 1 (q 4 I M + A T ) 

= ^ q 4'*'M ” ^ ^ q 4^M + ^ 
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if A is skew symmetric {which has been assumed) , and if the determinant 
of (q 4 I M ~A) is different from zero, i e , |q 4 I M ~A| 0 

From (3-2) , the characteristic equation of A can be determined, 

A 3 + w 2 X = 0 

Because a matrix satisfies its own characteristic equation , 

A 3 + w 2 A = 0, 
or 

a 3 

A = - — j , w f 0 (3-6) 

From (3-5) , again, 

q 4 X M + A “ M(q 4 T M " A) ' (3-7) 

or, using (3-6) and (3-4) , the left-hand side of (3-7) is 

a 3 

q 4 X M + A = q 4 X M ' “7 
w 

= ±2 (v 2 I m + 2q 4 A + 2A 2 ) (q 4 I M - A) , V f 0 

Comparing the above with (3-7), it is found that M can always be 
defined by 

M — (V 2 I M + 2q 4 A + 2A 2 ), v > 0, (3-8) 
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although it has yet to be shown that M is necessarily orthogonal when 
|q 4 I M -A| = 0. 

From the characteristic equation the eiaenvalues of A are found 
to be 0 , iw , and -iw For the case |q 4 I M ~A| - 0, must be an eigen- 
value of A, but, since q 4 is real, it must be identically zero In 
this case, (3-8) becomes 


M = 



w > 0, 


and, using (3-6), M can be solved for. 


M *= 



A=0 


A 2 =-w 2 


M 


(3-9) 


T —1 

In either case, M = M , and M is orthogonal This means that while 
(3-5) holds only for |q 4 I m - a| f 0, equation (3-8) is valid for all 
v > 0 


At this point it is necessary to examine the matrix M m more 
detail Since M is a rotation matrix, one of its eigenvalues (X^, for 
example) must be equal to unity The corresponding eigenvector coincides 
with the axis of rotation, since it is not altered m the rotation 
It is also known (Hildebrand, 1965, p 51) that 
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X 1 X 2 X 3 = l M l 

= 1, (3-10) 

and that 

\ ± + X 2 + X 3 = trace (M) * (3-11) 

A second matrix, M 1 , can be formed to perform the same rotation 
(Goldstein, 1950, p 123), 

cos ip sin tp 

M' - -sm ip cos ip 

0 0 

where ip is defined as the angle of rotation Since only a similarity 
transformation is involved (Hildebrand, 1965, p 54), 

trace (M) = trace (M 1 ) 

= 1 + 2 cos ip (3-12) 

From (3-9), (3-10), (3-11) and (3-12), 

2 COS ip = *2 + y- 

Solving for \ 2 and X 3 , and using the results of (3-9), 

A x = 1 , X 2 = e 1,fJ , X 3 = (3-13) 
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and, interestingly, it can be shown that the first eigenvector of M is 
simply 


S.1 



q 3 


From (3-8) , 


q 4 +q l" q 2 _q 3 2 ^ q l q 2 +q 3 q 4 ) 



2 (q 1 q 2 -q 3 q4) 


2.222 

q 4 +q 2 _q l’ q 3 


2 (q l g 3 +q 2 q 4 } 2 (q 2 q 3 _q l q 4 ) 


so that 


trace (M) = 


,2 2 

4q A - v 


Therefore, using equation (3-12), 


2qJ - 1 

COS iJj = 2 

V 


cos | = + — 


1 || 2 , 2 , 1/2 
sin f = +(1 - q 4 /v ) ' , 


2(q 1 q 3 -q 2 q 4 ) 

2(q 2 q 3 +q 1 q 4 ) 

2 2 2 2 
q 4 +q 3" q l- q 2 J ' 

(3-14) 


(3-15) 


where the signs in the last two expressions depend on the quadrant, 
which is defined by the first equation. 
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The matrix M can also be expressed in terms of the Euler angles 


cos to cos ft 
-cos I sin ft sin w 


cos ai sm ft sin w s xn I 

tcos I cos ft sm to 


M = 


-sin (c cos ft 
-cos I sm ft cos a) 


-sm w sm ft cos uj sm X 

+cos I cos ft cos 0) 


sm I sin ft 


-sm X cos ft 


cos I 


(3-16) 


A comparison of (3-16) with (3-14) yields the values for the q’s m 
terms of the Euler angles. 


< X i »ft— ii) \ 

q-L = v sin(j) cos (-j-) 

,1, /ft— UK 

q 2 = v sm(y) sin (-3“) 

q 3 = v cos ( 2 ) sin (-3—) 

,It ,ft+C0v 

q 4 = V COS(j) COS (-g— ) 


(3-17) 


Finally, it should be noted that the parameters p, a and t, used 
by Courant and Hilbert, and which will be used m section 3 3, are 
related to the Euler angles by 


a ~ 


ft-ai 7T 
2 ~ 2 


P * 


ft+tt 

~2 


T 



(3-18) 
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or, alternatively 


ft = a + p + ~ 

ir 

w = p - a - j 


X « 2t . 


3 2 The Transformation m Terms of a Set of Complex Polynomials , 


m,k 

b 2n 


<v 


Two new complex variables, and , are now introduced, and the 
coordinates x,y,z are mapped onto them by 


= x + iy 

2 

w 2 = -x+iy 


w 1 w 2 = z , (3-19) 

2 2 2 

where x +y +z =0 A transformation matrix, B, can, m turn, be con- 
structed to linearly transform w^ and w^, just as M transforms x, y 
and z. 


— — 


- — 

w i 

= B 

w i 

W 2 

i 

W 2 

_ 


_ _ 


( 3 - 20 ) 


The matrix B has two special properties (Goldstein, 1950, p 110), 


B B T = and |b| = +1, 
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such that if 


B = 




f 


then 

b l b 4 " b 2 b 3 = 1 
b l 5 l + b 2 b 2 = 1 
b 3 b l + b 4 b 2 = 0 

From these, and using (3-1), (3-14), (3-19) and (3-20), it can be 
verified that 


q 4 -lq 3 
. q 2 +lq l 


-g 2 +iqi 

q 4 +iq 3 


( 3 - 21 ) 


The new variables are - also introduced into the left-hand side of 
equation (A-ll) (through (3-19) ) and expanded m a finite binomial 
series , 


(-i)"(w 2 +w 1 t ) " = r (f)wf-*'(w 1 t) : 


2n 


= t n r ( f)w^y t ^ n 


Matching the coefficients of similar powers of t n with those in 
(A-ll) yields 
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(3-22) 


111 : H (r) - < 2 'L>"? + 'V'" 

n,m — n+m l z 


Making use of (3-20) and (3-21) , an expression analogous to (3-22) can 
be written as 


(i, v»® ■ e„>»rs n - m 


■ ( n«> (■/) (-q 2 -nq 1 -Kq 4 -iq 3 )s] 


n+m 


x [g4 + iq3 + (q 2 +i q 1 ) s ] n m 


(3-23) 


where s = w^/w^ 

-i v 

The next step is to introduce a new polynomial, (q^) , 

and a dummy variable, t, through two expressions for a particular 
generating function, G 2n (q v ,s,t), 


*2n 


(%■ = .*> - f f (f)v 2 »s3:”- k -"(q„)t3s k 

v D =0 k=0 3 


= [iq 3 (1-st) + 2 ^ (s+t) +q 2 (s-t) tq^^ (1+st) ] 2n , 

(3-24) 

where, although it is not proved here (Courant and Hilbert, 1966, p 542), 


4 


I 

v — 1 



m,k 

S 2n 


<V 


0 


The right-hand side of (3-24) can be expanded m a binomial series, and 
the coefficients of powers of t 3 equated with those on the left-hand 
side. 
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I v2ns 2n n ' k = [q 4 +iq 3 +s(q 2 +iq 1 ,l2n 3 C_ ^2 +;Lq l +S 1 3 

(3-25) 

Since s = w ]/ w 2 ' and letting j - n+m, the left-hand side of (3-25) is 


2n 

I v 

k=o 


2n Q m,k-n 
^2n 


(%) W 1 W 2 


The right-hand side can be identified with part of (3-23) , and, 
therefore, is equal to 


(1) (w 2 /v) 


-2n 


< 2n ) 

'n+m' 


Hence , 


( 1 ) 


Vm ( ^ " X S 2A k ^ n(q v )w l w 


2n-k 

2 


(3-26) 


From (3-22) , 


w k w 2n-k 

12 


( 1 ) 


H n,k-n ^ 

d n ) 


and inserting this into (3-26), followed by a change of summation 
index, gives 


( 1 ), 


n,m 


(R) 


n 

= I 

k— - 


( 2n ) 
v n+m' 

(2n » 

l k+n J 


S?' k (q ) 
2n KH v' 


( 1 ) 


H n,k<i> 


(3-27) 


The transformation formula for the spherical harmonics themselves 
can be easily obtained from (3-27), using (A-12) and the fact that R=r 
under pure rotation The transformation from the old harmonics to the 
new ones, in terms of the S? ,k, s, is 

Vn r 
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P™(cos aje 111 ^ 


n 

I (-D 

k*-n 


k-m (n-k) 1 


(n-m) ! fa 2n 


• m ' k P k (cos 9)e lk * 


(3-28) 


in k 

3 3 Evaluation of the Polynomial as a Function of the 

Rotation Parameters 

In section 3*2, the transformation of the spherical harmonics was 
formulated m terms of a set of polynomials, (q^) ' s , which have 

not, as yet, been evaluated in this section, they will be determined 
as functions of the Euler angles associated with the rotation of 
coordinates 

Introduction of p , a and t, as defined by (3-18), allows the g's 
to be expressed as 


q^ =-v sm o sm t 

^2 “ v cos 0 sin T 
= v sm p cos t 

q 4 = v cos p cos t 

The generating function G 2 n ^ v ;s,t ^ can then be written (from the 
last part of (3-24) ) as 


G 2n (p,a, T/ S,t) = v 2n [cost e lp +s sint e 


-1 a 


-t sinr e ia +st cost e lp 3 2n 


This can be simplified by letting 
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-i (p+a) 


-l (p-a) 


(3-29) 


s* = se 


t* = te 


so that 

G^ n - v 2 n e l2pn [cost (l+s*t*)+sxnx (s*-t*}] 2n 

But, from the first part of (3-24) , the same quantity can be expressed 
as 


G 


2n 


n n 


I I 

m=-n k--n 


2n x 2n r ,m,k, n+m k+n 
n+m v s 2n * s 


n n 

l 1 

m=-n k=-n 


,2n , 2n l (k+m) pi (k-m)a . *n+m c *n+k 2ipn_m,k 
n+m 2n 


Equating these two expressions yields 


[(l+s*t*)cos t + (s*-t*)sm x] 2n 


= i i 

m=-n k— -n 


.2n . m,k 1 (k+m) p 1 
m+n' b 2n e e 


(k-m) a t *n+m g *n+k 


(3-30) 


The left-hand side of (3-30) can be expanded in the familiar 
binomial series, and is 


l t*P +n ( l 2 ^ ) (s* cos T-sm T)P +n (cos t+s* sm t) n p 
p=-n p 

The coefficients of powers of t* n+ ™ xn the above series and in (3-30) 
are equated and then multiplied by cos 11 m r sin m+n i: , giving 
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cos n " m T sm n+Itl T I S *n+k s m,k e i(k+m)p el (k-m)a 
k=-n 2n 

* 

/ 2 , . ,n-m, * , , 2 .n+m 

— (cos t + s*cos t sm t) (s*cos t sm -1 + cos t) 

(3-31) 

To further simplify equation (3-31) and what follows, let 
2 

a = cos x 
b = cos t sin t 
X = s *b , 

and apply Maclaunn's series in powers of x to the right-hand side 
of (3-31), 

l j£ ^ (3-321 

p*° 1 P a * P J x ,o 

The evaluation of the quantity within the bracket is accomplished 
using a result of Hobson (Hobson, 1965, p 125), 

~k ~k 

— <f> (a+x) = — v 

_ 9x k j x _ 0 a** 

In this case, let 

4> (a+x) = (a+x) n-m (a+x-l) n+m , 
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so that (3-32) becomes 


2n 

l £ 

p=0 P 


*P 


tos^T sm p T 


d ( cos^t)^ 


[cos 2(n - m) r (cos 2 T-l) n+m ] f 


since 


§_ = l_ 

da 9a 


Substituting the above series into the right-hand side of equation 
(3-31) results m a power series in s* on both sides of the equation, 
so that the coefficients of s* n+k can be equated With this, 
can finally be solved for explicitly. 


, k > \ 

S 2n 


-x u 


r 2 (n-mL , 2 ix n+m, 

[cos (cos t-1) ] 

(3-33) 

The derivatives m (3-33) can be related to the Jacobi polynomial 
J n , which, in turn, can be expressed as a hyper geometric series, 2 F i 
In general (Courant and Hilbert, Vol 1, pp 90,91), 


J n (p ' q ' x) H Yr 3L [x« +n “ 1 (l-x)P +n - <I ] 

q (q+l) . (q+n-1) - n 


= 2 P 1 (-n,p+n,q,x) q > 0, |x| < 1 

Making the following identification of indices. 


n+k n, n+l-k -*■ p+n, 1-m-k + q. 
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equation (3-33) becomes, for m+k £ 0 and t > 0, 


S^ k = ( _ 1) n+m ( n-m )e -i(m+k)p e -i(k-m)( rcos -m-k T gin m-k T 


x 2 f 1 (-n-k,n+l-k, l-m-k;cos x) , 


(3-34) 


where 


F (d g y ?] = i + x + - ( -a + 1) . B .^ +1) x 2 -*- 

^(“/BfY/XJ _ 1 + jry x + i 'yCy+l) x 


(3-35) 


It can be shown that 


= m,k n .m+k c -m-k 
s 2n " l "- u fa 2n 


(3-361 


so that for m+k 0 and x > 0, 

m,k = 1) n+k ( n+m )e -i(m+k)p e -i(k-m)a cos m+k T sin k-n» T 
zxi n-K 

x (-n+k,jvfl+k, 1+m+k, cos 2 t) (3-37) 

In the special case when 1=0, can be evaluated directly 

from (3-31) , all the polynomials are zero except for k = m, where 

S^ (T = 0) = e“ l2mp 

m k 

Using (3-18) again, the expressions for are written m 

terms of the Euler angles ft, I and w, 
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<,m,k 

s 2n 



(„ 1 ) n+ ra ( n-in) e -i(m-k ) Tr/ 2 e -i(ii«o+kn ) cos -m-k I sir m-k I 

x 2F 1 (~n-k,n+l-k,l-m-k,cos 2 i) , m+k £ 0, I f 5 0 , 

( _ 1) n+k ( n+m )e -i{m-k)Tr/2 e -i(mtfl+kJ2) oos m+k I sin k-m I 
n— k * ^ 

x 2 F 1 (-n+kyn+l+kfl+m+kjcos 2 j) , m+k £ 0, I ^ 0 , 


(3-38) 


(3-39) 


-im(w+n) 

0 


k = m, 1=0, 



k ^ m, 1=0 (3-40) 


Ttt V 

when k assumes the value of zero, S 2 ^ is independent of ft (and 
in k 

vice-versa) , and S 2 ^ reduces to a function involving a Legendre func- 
tion Prom (3-39) for k»0 and m > 0, 


c™/° 

b 2n 


t , n -iti/n+iru “imw 
(-1) l ( )e 


cot m 2 ^^ (-HfH+1 , 1+m, cos^ 5*) 


(3-41) 


But m this case the hyper geometric function is of the form (Abramowitz , 
1965, p 332) 


F(-n,n+l,l-p,Li) = r ( 1 -p) (1 - P //l P%(x) 
2 (l+x) p/2 n 


Letting x = -cos I and m = -p, this becomes 


F (-n,n+l, 1+m, cos 2 |) = m’ tan m | P n ra (-cos I) 
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Then, using (2-4) and (2-5), this is 


F(-n,n-M,l+ra,cos 2 j) = (-l) n m l tan m j "7 ^4-mj ’ p ^ cos 

(3-42) 

Substituting (3-42) into (3-41) yields the special result for k - 0, 
which is valid for all m, 

S 2 ^ - i n i - e P n (cos I) (3-43) 

Due to the lack of standardization of the notation m the litera- 
ture, Table 3-1 has been included to assist anyone researching the 
subject 
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Table 3-1 Notation Used by Several Authors to Represent the Same 
Quantities 



Lee 

Courant, 

Hilbert 

Aardoom 


Jeffreys 

Gold- 

stem 

Subscripts 


order 

LL. 

i 

i 

i 

n 

i 

i 

* n 

i 

i 

i 

i 

_L. 

a 

i 

i 

i 

L- 

i 

i 

i 

i 

degree 

r 

m 

1 

1 

l 

— j 

j s 

1 

I 

m 

1 

1 

1 

i 

i 

summation 

k 

I 

1 

1 

j 

r 

i k,* 
1 

r- 

i 

I 

r 

s 

f-’ 

1 

i 

l 

— j 

i 

i 

i 

Polynomi als 


I 

1 

1 


r 

i 

| . 

i 

i 


1 

1 

I 

l 

i 

i 

i 

t 

associated 

Legendre 

P m 

n 

1 

1 

1 

1 

1 


1 

! p s 
i n 
f 

i 

i 

i 

_i 

(*-m) 1 

pm i 

rP Jl ! 

i 

i 

i 

i 

transfor- 

mation 

r 

qltt t 

s 2n i 

i 

i 

t 

■& r 

1 

| ^ 1} s+Ag n+s,n+& 
2n 

i 

i 

i 

i 

i 

i 

i 


r 

i 

i 

i 

i 

i 

f 

i 

i 

i 

i 

i 

Euler Angles 


1 

1 

1 


i 

i 

i 

i 

i 


i 

i 

i 

i 

J 

* 

i 

♦ 

longi- 
tude of 
ascending 
node 

fl 

J 

l 

l 

l 

l 

.1- 

a +P+ | 

i 

i 

i 

i 

j 

i 

i 

! 

n, n+2 

i 

i 

i 

i 

i 

i_ _ 

! 

+ ! 
i 


T 

i 

I 


r 

i 

i 

1 

l 

1 


i 

i 

i 

~! 

i 

i 

inclina- 

tion 

I 

l 

l 

l 

l 

2t 

j e 

i 

1 

! 

1 

1 

I 

i 

i 

i 

i 

i 

e j 

i 

argument of 
perihelion 

w 

i 

l 

i 

l 

i 

l 

TT 

1 V — — 

j X 2 

l 

I 

“1 — 
1 
1 
1 
1 
l 
1 

~ TT 

“'*-2 

■" i — 

i 
i 
i 
i 
i 
i 

* I 

i 

i 

i 

Old Coordinates 


l 

i 

[ 

i 


1 

i 

I 

• 

1 

1 

1 

— «_ 


i 

r 

\ 

i 

i 

i 

co latitude 

0 

l 

l 

i 

4* 

1 TT A 

j 

l 

l 

1 

0 

— f— 

i 

i 

i 

i 

i 

i 

i 

longitude 

$ 

i 

I 

l 

X 

! x 

i 

-r-- 

i 

i 

X 

■ — ** 

1 

l 

j 

i 

radial 

distance 

r 

*r 

i 

i 

i 

i 

i 


i — 

! r 
i 
i 
i 

i 

i 

i 

i 

i 

i 

r 



1 

l 

l 

l 

l 

( 

i 

i 

i 

i 

i 

i 

New Coordinates 


i 

i 

i 


i 

i 

i 

i 

i 

i 


1 

I 

I 

i 

i 

i 

co latitude 

J 

a 

■f 

i 

i 

r 

♦ ’ 

! l-r 

' 2 9 __ 

i 

i 

i 

0 1 

1 

I 

1 

1 

I 

* 

i 

i 

longitude 

B 

( 

! 

x r 

T 

! x' 

_L _ — 

~j 

i 

i 

_t 

X 1 

| — 
1 
1 

7 

i 

1 

radial 

distance 

R 

1 

1 

1 

1 


i r 1 

i 

i 

"r- 

1 

1 

1 

I 

r 1 

— 1— 

1 

1 

1 

1 

— — -j. 

[ 

1 

i 

t 
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CHAPTER 4 

TRANSFORMATION OF A SERIES EXPANSION 
IN SOLID SPHERICAL HARMONICS 


4*1 Explanation of the Transformation Formulas 

In chapters two and three , transformation formulas were developed 
for transforming solid spherical harmonics under changes of coordinate 
axes Two classes of harmonics were treated because a general expansion 
of a harmonic function will involve either, or both, kinds A func- 
tion, V, for example, m the old coordinate system, is of the form 

co s 

V(r,e,4>) = l l P*(cos e) j (;) s+1 [c cos t<i>+S sm t<i>] 
a s=0 t=0 ' ^ 

S J 

+ (|) [E st cos t4> + F gt sin t<j)]j, (4-1) 

where G, M and a are all constants, along with the coefficients C 
S , E st and F st ^ see table As ^plied by the notation, V can 

represent a potential function m free space (Kellogg, 1953, p 218) 

If V is a gravitational potential, G is the universal gravitational 
constant and M is the mass of the gravitating body The form of the 
expansion m (4-1) is probably the more familiar, but a slightly dif- 
ferent, more symmetrical form will be used here, 

00 s a- i S+l i , 

V(r,0,<f>) = f£ l l JpJ(cos 0)| (5) 

+(|) S Ee[(T st+1 W st )a- lt *]] , (4-2) 
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where Re stands for the real part of quantity enclosed by brackets 


Old Coefficients 


New Coefficients 



Hence , the coefficients m the two representations are related simply 
by 



t > 0 



(4-3 ) 


Analogous relationships hold for the other four sets of coefficients m 
Table 4-1 
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After translation and/or rotation, the new series will be, m 


general. 


V(E , =,,♦). p i t »>=*«, irww- - *i 


n=0 m--n 


+ (!l nR «! l0 »» +5 »m«' V »» t5 »m>) e " :i ” e !|- ' 4 - 4 > 


Once the coefficients are determined, they may be converted to the asym- 
metric series coefficients by (4-3) , if so desired 


Asymmetric Form Symmetric Form 



Nothing has been said, thus far, regarding the convergence of the 
various series expansions The original series, whether it is xn 
powers of r 11 or r n or both, is assumed to converge For special 
cases, the constant "a" can be associated with a characteristic length 
of the gravitating body, and convergence of the series linked to the 
ratio of r/a 

The primary concern here is the convergence of the resultant 
series In the case of a straight rotation of axes, as was 
seen in 3 2, there is no problem, as the A gtl B st l s transform into the 

J nm' K nm' s and fche T st' W st' s into the U rnn' V nm’ s Wlth a translation, 

however, the T . ,W . f s transform into the U ' s, while the A . ,B . T s 

st st n m nm st st 
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transform into U .V T s if R<p, or J / K 1 s if R>p (See figure 4-1 
nm nm nm niu 

and table 4-2) 

In other words, the exact form of the series m the latter two 
cases depends on the ratio of the lengths of both the field and trans- 
lation vectors The critical sphere separating the two regions of con- 
vergence is of radius p, centered on the new origin, and passes through 
the old origin 



Figure 4-1 Critical Sphere for a Transformation Due to Translation 

-s-1 

of a Series Involving Powers of r 

The various cases are treated separately m the following four sections 
of this chapter 
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- s 

4 2 Translation Transformation of a Series Involving Powers of r 
to One with Powers of R n 

In this section it is assumed that a function, V(r,0,4>), has been 
expressed m terms of convergent series of solid spherical harmonics 
of the following form 

v = l f (|)V(cos 0)Re[(T .+iW .)e" xt *] , (4-5) 

a s =o t=-e a s st st 

where the T st ' s and W st 's are known With a pure translation of the 
origin the same function is expressed in a new series m the new 
coordinate system. 


V (R,a, B) 


GM y 
a L 


R ^ 

- I <i> 

n=0 m=-n 


pJ(cos 


(4-6) 


where the U 1 s and V _ 1 s need to be evaluated 
nm nm 

Equation (4-5) can be written as 


V = Re 




s . 

l Pg (COS 

t=-s s 


e)(f) s (T st +iw st )e 


-it<J> 


Using (2-16) and (2-17) , which transform the harmonics to the new 
coordinate system, this becomes 


V = Re 


GM 


i i <*. t +i "*)* 

s=0 t=-s 


. s l k , 

-s r r|-k+s+t 

L Z I L- 

k=o ^| k ; t _ s 


x ( k+A )pS_ke_ 1 (t " &)Yp s-k (cos X5Rke ^P^Ccos a) | 

(4-7) 

Switching the inner two summations to the outside, and changing nota- 
tion, (4-7) can be written as 
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co n co s-n+m 

V = |^ i i i i O a ' n (|) s “ n « n 
a n-0 m=-n s=n t=n+m-s a 


X P™(cos ct)Pg"^(cos X)Re[(T st +iW st )e -l{t m )Y e -imP] 

(4-8) 




The real part of the bracketed portion of (4-8) is 


cos mg[T st cos (t-m)y + W gt sm (t-m)y] 

+sm mpt-^^sin (t-m)y + W gt cos (t-m)y] (4-9) 

Inspection of equation (4-6) , the other expression for the poten- 
tial, shows that the harmonics of degree n, order m, in the same 
variables, must match those m (4-8) The coefficients are then 
equated, where the real portion of the brackets in (4-6) is 


U cos mp + V sm m3 
nm nm 


Thus, (4-8) is 


U nm 
V nm 

T st cos (t-m)y + W gt sin (t-m) y 
_-T st sin (t-m) y + W gt cos (t-m) y 


» s-n+m 

■ I l 0<i> 

s=n t=n+m-s n+ITl a 


s-n 


P^(cos X) 


or using (2-4) and regrouping the quantities within the brackets. 
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oo 


U +iV 
run nm 


1 

s=n 


s-n+m 

I 


t=n+m-s 


(-D^ts+t) 1 ,£>s-n 

(m-n+ s - 1 ) r (n+m )’ 


X p“~* (cos X)e l(ltl “ t)Y (T st +iW st ) (4-10) 


4 3 Translation Transformation of a Series Involving Powers of 

r s 1 to One with Powers of R n , R<p 

In this case, the coefficients A and B are assumed to be 

st: st 

known m the expansion 


V(r,6 ,4>) 


GM y 
a t 


! (|) S+ 1 pJ(cos 
s=0 t=-s r 5 


e)Re[(A st+ 1 B st )e- 1 ^] 

(4-11) 


The translation of the origin transforms the series into the form of 
* * 

(4-6), with the U *s and V 's unknown 
nm nm 

The method of solving for the coefficients is analogous to that 
of section 42, using equations (2-22) and (2-23) to substitute for the 
"old" solid harmonics m (4-11) The result is 


* * 

U +iV 
nm nm 


V f (_l) n+m+t (n+s-m+t) ' 
s=0 t=-s (s-t)Mn+m) ' 


X 


( 4 - 12 ) 


4 4 Translation Transformation of a Series Involving Powers of 

r to One with Power of R n R>p 

Here the coefficients m the old series (4-11) are known, but 

the J 1 s and K ' s are to be determined for the new series 
nm nm 
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V(R,oc,|3) = l l (f) n+ Vj(cos alSeUJ^+iK^le" 1 * * * ” 8 ] 
n-0 m=-n 

(4-13) 


Again , the procedure is the same used in 4 2 and 4 3, equations 
(3-22) and (2-23) provide the transformation formulas for the solid 
harmonics The result is 


J 


+1K 

nm nm 


r yjn-fm-s (-l) n S (n-m) 1 
t n JU-s (s-t) (n+m-s-t) ! 
"° j s+m-n 


* ( i >n Sp n-l (cos X)el(m " t)Y(A st +iB st ) (4-14) 


4 5 Rotation 

A slightly different approach is taken m the case of a rotation, 
since the new harmonics have been determined in terms of the old ones 
(3-28) , instead of vice-versa, as for the translation The two 
approaches are equally correct, the appropriate choice depends entirely 
on the direction m which the harmonics are transformed 

Because the magnitude of r remains unchanged under rotation, there 
is only one form for the transformation, as opposed to the three for the 
translation Series involving A st r B st 'z and T st ' W st ls are * therefore, 
transformed by the same formula The first case is worked out here 
The A^'s and B st ' s are assumed to be known in the expansion 


v(r,e,4>) = I 


l (f) S+1 pJ(cos 6)Re[{A st +iB st )e 

3-0 t “" S (4-15) 


But, as pointed out in 4 1, V can be considered as a potential func- 

tion Hence, it can be represented by an integral over the body 

generating the potential. 
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V (r ,0 , <j>) - G f dm ' - , 

M |r~r T I 

where the primes indicate variables of integration, dm r being an element 
of mass, m the case of a gravitational potential 



Figure 4-2 Definition of the Variables of Integration, r r ,9 , ,<i> , ,tj>', 
for an Element of Mass at Q and a Field Point at P 

Using the law of cosines m terms of the angle between the field 
vector and the integration vector, 

Ir-r 1 | = (r 2 +r 1 2 -2rr 1 cos 

Then, applying the generating function for Legendre polynomials (2-9) , 
followed by the addition theorem (2-10) , and using (2-4) , the potential 
is 
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v(r,6,4>) 


6 



-n 


(n-k) 1 
(n+kT r 


pj^ (cos 


0) 


j< P^(cos , 


where 



r 


< X 



> 1 


(4-16) 


(4-17) 


Equating the coefficients of each harmonic P^(cos 0) m (4-15) and 

s 

(4-16) yields 


a 

M 


J dm’r'pJ(cos 
J H s s 


a l \ A Xfc<j> ' _ (s+t) 1 
0 “ (S-t)"' 




(4-18) 


After the rotation/ the "new" potential/ V(R/a,3)r 13 of the form 
of (4-13) But it, too, can be written m terms of a volume integral 
m exactly the same way as V(r,e,(i>) Thus, 


V(R,a f 8) 


/ dm 1 


00 n 

1 R A l 

n=0 m=-n 


(n-m) 1 
(n+mp 


P™(cos a') 


x (cos a) Re [e~ lin ^ ® ^], 


or, equating the coefficients of the harmonics, P™(cos a) , 


J +1K « S 

nm nm M 


{hsK ! “ s 
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from (3-28) , this becomes 


Substituting for P™(cos a')e lln ^ 


J ™ WK ™ - s X,, 1 ' 11 ’'"”' tSht- s ?A k Jj^'K ^<“ s 6,) 


ik(Jj * 


(4-19) 


But , since R=r and R f =r f , the integral (4-18) can be substituted for 
the integral in (4-19) , resulting in 


, +lK ™- l {tar s ?A k(A „k +u W 


nm nin 


k--n 


(4-20) 


As pointed out earlier , the relationship between the U nm , v nm ' s 
and the 1S same / 


U +iV = T (-1) 
nm nm L v ' 


k=-n 


k-m (n+k) 1 0 m f k /r71 * 

Tn+mf' S 2n (T nk +lW nk ) ' 


(4-21) 


and note that only coefficients of the same degree (n) are used to 
determine the new ones 

Finally/ as a matter of completeness, it is possible to indirectly 
solve for the old spherical harmonics in terms of the new ones (under 
rotation) using (4-20) and (3-28) In chapter three the harmonics were 
transformed in the opposite order - the new ones m terms of the old 

With the method developed in 4 2, it can be easily verified that 


Pg (cos 


G)e 


it<j> __ 


s 

= i (-D 

m=-s 


t-m 


(s+t) 1 
(s+iu) 1 


S 2 s ,t? s(- S 


a) e 


im6 


(4-22) 


Comparison of this with (3-28) reveals their similarity, except that in 
this case the complex conjugate of is involved, and it is summed 

over the first superscript as opposed to the second m (3-28) 
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CHAPTER 5 

THE FORCE BETWEEN TWO HOMOGENEOUS HEMISPHERES 


5 1 Potential Due to a Hemisphere 

One of the primary objectives of this thesis is to demonstrate 
the applicability of the transformations developed in chapters 2-4 to 
problems involving volume integrals of bodies possessing rotational 
symmetry As an example, the force between two bodies A and B can be 
calculated from 


F = 


p d 


/ 

volume 
of B 


vv ax. 


where is the mass density of A, and V is the gravitational potential 
of A If B is homogeneous this can be expressed, using the gradient 

theorem, as a surface integral involving V, 


F 


p d 


/ 

surface 
of B 


VndS, 


(5-1) 


where n is the outward normal to the surface 

In the evaluation of the integral m (5-1) it is desirable to 
make the integrand as simple as possible This is accomplished by se- 
lecting the appropriate type of coordinates - m this case spherical - 
and orienting the coordinate system to take full advantage of the 
symmetry of body B But V also has to be determined m the coordinate 
system best suited for body A Hence, once V has been calculated, the 
coordinates are translated and rotated for the evaluation of the inte- 
gral in (5-1) As shown m chapter 4, the potential V is expressed m 
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the new system with a set of transformed coefficients ♦ 

This method is ideally suited for two homogeneous hemispheres, due 
to their simple shapes The potential of a hemisphere along its axis 
of rotation can be found from (see Figure 5-1) , 


v = p a G 


f f* f a dr 

'*'♦-0 * 4 » ir /2 * R =0 ** 


(5-2) 



I 

I 


Figure 5-1 Calculation of the Potential Due to a Hemisphere 
at a Point on Its Polar Axis 


Evaluation of (5-2) yields 

r 

V(z) = | irp a G d 


1 + f 


+ ( l )3 


[1 + (^-) 1 

ct 


2 3/2 


(5-3) 
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The off-axis potential is determined by expanding (5-3) m powers of 
(z/a) n or (z/a) n 1 and substituting (r/a) n P n (cos 6) or (r/a)~ n ~^P n 
(cos 6), respectively (MacMillan, 1958, pp 360-362) This gives 


V(r,cos 0) = |rrp d Ga 2 | | - | | P i( cos e > + ( I )Zp 2 (cos e) 


r\2 T 


■| ( i )3p 3 (cos e > + 3 l 


(-l) n (2n-3) ' 
n=2 2 2n_1 (n+l) ' (n-2) ' 


r 2n+1 ) 

x( S> P 2n + l (cos 0 > J' 


0 < r < a, (5-4) 


or 


2 a 3 ,a,2. 


V(r,cos 0) = jiTp^a j = - ~ (^^(cos 9) 


1 ,a, 4^ , , “ (~l) n (2n-l) • 


♦ 16 V«» 5 6 >- 3 l .in*!, — ,,, 

n=2 2 (rt+2) ' (n-1) ' 


v ,a. 2n+2„ , 

X( ? } P 2n+l (cos 0) |' 


r > a 


(5-5) 


5 2 Force Between Two Hemispheres 

The force exerted by a hemisphere on a second identical one is 
calculated using (5-1), (5-4) and (5-5) The center of the face of the 
second is located by the vector p relative to the first (Fiqure 5-2 (a)) 
The force vector F is determined relative to the new coordinate 
system £, n, C, or R, a, $ In this system n dS is 
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n dS 



0 

0 

-1 


R sm a dR d3 


sxn a cos 3 
sm a sxn 3 
cos a 


R sxn a da d3 


flat surface 
(a = tt/2) 


spherxcal surface 
(*=a) , 


so that (5-1) xs comprxsed of two xntegrals 


F = F f + F s 


/ 


P* J V (a = tt/2) 
a flat 
surface 


0 

0 

-R 


dR d3 


+p d a 


/ 


spherxcal 

surface 


V(R-a) 


sxn a cos 3 
2 Q 

sm a sxn 3 
sxn a cos a 


da d3 


(5-6) 


From (4-4) the potentxal of hemxsphere A xn the two cases , after coordx- 
nate transformations , xs 


co n +l 

V(ct = tt/2) = t l P™(0) j(|) [J^cos m3 + K sin m3] 

n=0 m=~n * ' L 1 1 


+ ( l> n[0 M , +5 ™>' 3 “ "8 + < V m Am ,sln "«! , (5-7) 
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and 


V (R=a) = — l l P m (cos a) [ ( J +U +U )cos rag 


n=0 in=-n 


n 


nm nm nm 


+ ^ K nm +V nm +V nm^ Sln »« 


Although the extreme limits of the integration variables in 
are specified. 


0 < R < a 


(5-8) 

(5-6) 


0 < a tt/2 
0 < 6 £ 2 ir, 

the particular limits corresponding to the region of convergence for 
determination of the coefficients m (5-7) and (5-8) is somewhat com- 
plicated, depending on the relative positions of the two bodies 

In general there are five distinct cases for evaluating (5-6) 

The first is the case when the two hemispheres are face-to-face with no 
separation, i e , p = 0^ Since all surface points relative to the old 
origin, r s , are equal to the radius, and since all R g are greater than 
p, the integration of (5-6) involves only the coefficients of the 

U ,V lc * or the J__ ,K 1 s , but not both 
nm r nm nm nm 

The second case is distinct because for every surface point, either 

r g < a or R s > P The integral, therefore, involves only u nm ' V nm ,s 

and J ,K 1 s This includes all p up to P-, , where p- is the radius 
nm nm l i 

of the sphere centered at the new origin, passing through the old origin, 
and reaching out to the first intersection of the surface of hemisphere 
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Either Flat or Spherical Sur- 
Valid for face Integrals, or Both, In- 

Case Criteria These P volve These Coefficients 


I r s =a and p=0 
R s >p 


U ,V 1 s or J ,K 's 
nm 7 nm nm 7 nm 


II 


r <a or 
s 

R >p 
s y 


or 


U ,V * s and J -K T s 
nm 7 nm nm 7 nm 


III 


general 

case 


P 1 < P<P 2 


U ™' v ' s J ,K 's and U ,V 1 s 
nm nm nm 7 n m nm nm 


IV r g <a or P 2 <P<P 3 

V p 


U ,V 's and U .V 's 
nm nm nm' nm 


V all 


P 3 <P 


* * 

U ,V ' s 
nm 7 nm 


r >a 
s 


Critical p T s 


a/2 < p 1 < a/ /2 


a £ p 2 £ 2a 


a < p 3 < 2a 


Table 5-1 Summary of the Five Cases Involved m Evaluation of the 
Force Integrals 
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B with the sphere of radius a centered at the old origin (Figure 
5-2 (b) ) Examination of the geometry shows that 0 < p <_ and 

a/2 < < a//2 

The third case will be discussed after cases four and five The 

fourth case is special because it includes integrals of only the 
* * 

U ,V 1 s and U .V *s From Figure 5-2 (c) it can be seen that 
nm nm nin nm 

hemisphere B must lie within the sphere of radius a centered at the 
old origin, or within the sphere of radius p centered at the new ori- 
gin In this case either r g < a or R g < p* From Figure 5-4, 

P 2 < p £ P 3 , where 

a £ p 2 £ 2 a and a <_ <_ 2 a 

Due to the overlapping of the possible limits of integration, this 
case may not always be applicable 

The fifth case applies when all R g £ The integral, therefore, 

'k k 

involves only the coefficients U .V 's, and p 0 < p 

nm nm 3 

Case three is the general case and includes everything not covered 
m cases 1,2,4 and 5 In order to evaluate integrals over both the flat 
and the spherical surfaces , all the coefficients are used 

5*3 Evaluation of the Force when the Limits of Integration Have No 
Azimuthal Dependence 

For a particular orientation of the hemispheres the force can be 
calculated after determining the appropriate limits of integration from 
the geometry In this section it is assumed that there is no dependence 
on 3 / the azimuthal angle This includes two cases which will be worked 
out m detail - separation of the hemispheres along the z axis , without 


64 



A /Ji 

rotation, and a case involving rotation when £ 2ai z , where £ z repre- 
sents a unit vector m the positive z direction In this case, $ can 
be integrated from 0 to 2 tt , and (5-6) becomes 


u = 


f f 2v GM 
1 ''r*'s=o a 


R 


co n 

i 1 

n-0 m^-n 


pj(0) 


( 2 ) 

r 


n+1 


[J cos m3+K sm m3] 
nm run 


+ (f) n [(U +U )cos mB+(V+V)sin m3] 
a nm nm nm nm 


i ? dp dR. 


All of the terms involving m3 integrate to zero, except for in^O, leaving 


✓v /*R n p 

F- = -2irp,GM l V (0) / maX [ (|) J +(-) 

— f *a — C n •/ ^ L R nm a 


n+1 


n=0 1A R==R 


(U +U ) ] dR 
' nm nm 7 J 


mm 


(5-9) 


The polynomials P n (0) are obtained by Robin (1957, p 75), 


n-0 


p 2n (°) = 


(-l) n (2n-l) 
2 2n_1 n 1 (n-1) 


n=l ,2,3 


(5-10) 


where the polynomials of odd degree are all zero Those of even degree 
can be evaluated from 


P 2n< 0) 


(2n-l) - 
2n 2n-2 


( 0 ) , 


n > 2 


With this, equation (5-9) becomes 
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Ff = -2 1 rp d GMx ? {RJ 00 + £ (U 00+ U 00 ) 

R max 

R mn 
(5-11) 


00 

+ l * 

n=l 


2n 


( 0 ) 


(-) 


2n-l 


2 (n+1) 


2n0 


l-2n 


<i> 

2 (n+1) 


<WW 


Over the spherical surface (5-6) is 


F = a p . 


> 00 n 

J J JT I I F™(cos a) [J +U +U )cos m$ 

a*p n=0 m=-n n nm nm nm 


+ (K +V +V ) sin mBJ 
nm nm nm 


sm a cos 3 
2 

sm a sm 3 
sm a cos a 


dS da 


This time all of the terms with [m| >1 drop out when 3 is integrated 
from 0 to 2ir Consequently, 


F g - 2irp^GMa £ J* da 


rp n ( cos o)sin 2 a(J nl +U nl +U nl ) 
P n Ccos a ) sm 2 o (K n i + ^ nl + V nl ) 


P n (cos a) sin a cos a <J n0 +U n0 +V n0 > 


(5-12) 

When the hemispheres are separated along the z axis , or for rota- 
tions when p > 2a, equation (5-12) is integrated from 0 < a £ ir/2 
This involves two separate integrals of the form 
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X 


ir/2 


1 2 
P„ (cos a)sm a dot 


(5-13) 


and 


rf r/2 

/ P (cos a) sxn a cos a da 
•'0 n 


(5-14) 


The integral in (5-14) is evaluated by Robin (1957, p 26) If C n is 
defined as the integral (5-14) , then 


1/2 


2n = 0 


1/3 


2n = 1 


C 2n 


1/8 


2n = 2 


(5-15) 


- (~l) n (2n-3) ' 
2 2n_1 (n+1) ' (n-2) 


n “ 2,3 , 


or recursively. 


r = 3 ~ 2n r 
2n 2 (n+1) ( '2n-2 ' 


n > 3 


The integral in (5-13) can be expressed as 


r 


(l-v 2 ) 1/2 P^(v)dv, 


(5-16) 


but from (2-3) it can be seen that 
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- -(l-v 2 ) 1/2 e„<v). 


so that (5-16) is 


-Jo iv V v)dv + f 0 h V v)dv 

Evaluation of the integral above yields 

f\ l-V 2 ) 1/2 P 2ii (v)d» . P 2n (0) - 2C 2n 

The force due to the integration over the spherical surface 
(5-12) can now be evaluated, ana is equal to 


i w u«Ti + “ii> 


£ s = 27rp d GMa< 


- |« u « u 4 u ) 


^00 U 00 u 00 . lfj , n * , 

• 2 + ~ 2 ~ + ~2 + 3 (J 10 + U 10 + U 10 ' 


<P 2n (0) 2C 2n ) (J 2nl +U 2nl +U 2nl ) 


+ l ( p o„«>)-2C 2n )( K 2nl+ V 2nl 4 2nl ) 


(5-17) 


C 2n (J 2nO +U 2nO +U 2nO ) 


5 4 Calculation of the Force for Two Special Cases 

The results of 5 3 are applied to two cases - when the hemispheres 
along the polar axis, and when a rotation of hemisphere 
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B is performed in the Y-Z plane through an inclination X (Figure 5-2) 

For these positions case 4 (from Table 5-1) does not apply, and the 
first case is a derivative of the second, so that only one of cases 
two, three and five needs to be evaluated, depending on the value of 
p Because the rotation will be carried out for p = 2a, only case three 
is involved The equations developed m 5 3 are therefore applicable, 
as the integration over a and £ are independent of the translation and 
rotation 



(a) (b) 


Figure 5-3 Positions of the Hemispheres for the Two Special Cases 

Except m the case of the rotation, the force will have only a £ 
component , hence the total force may be broken down as 

F = (F + F + F,p 4* F- + F_ )i + F 
“ S 2 S 3 *2 t 3 £ 5 ~ S 5 

The geometry dictates that the critical values of p are p^ = a//2 
and p 2 = a, and the limits of integration for the various cases are 
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R = 0, (a^p 2 ) 1 ^ 2 and a Thus, from equation {5-17) and Table (5-1), 


F - ir Mp d a[J 00 + 3 J 10 + 2 T C 2n J 2nO ] ' 
^ n=i 


F = F 
s 3 S 2 


F = TrGMp d a 
5 


4 1 11 ! +2 i c ,J 


4 $ f, 2n 2n0 

3 V 11 ] n_1 

O0 

U nn + |L+ 2 I C 0 tL A 
00 3 10 *•- 2n 2n0 

n=l 


Similarly/ the forces evaluated over the flat surface are, from (5-11) 
and letting R^ = (a 2 ~p 2 ) 1//2 , 

V -*“P4‘ ho' 1 - 3 J 10 


R, 2n+2 

» P,„(0) * 2n-l P 5ri (0)(-i) 

" 2 £, J 2n0 [C 2n + - 2n-l"" 2n+2 U 2 

n—± 1 


F j - -* s »Pa a | -| J 10 + f 1 - 2 f> J oo + ( r ,2 °oo 


2 R, 2 


+ [ (|) t (jl ) Uqq+ 2 y 


R.. 2n+2 

P, n (°) . 2n-l P„(0)(-i) 

C 2n J 2nO Siv 1 ! J 2nO tl-( p ) ]+ 2n+2 U 2n0 


P 0 „ (0) ^ 2n+2 R. 2n+2 * 

+ r (£) _(JL) i5 

2n+2 u a ; l a ’ JU 2n0 


( 1 * ” * P 2n (0) ) 

r f- = - 27 rGMp d a j-| U 1Q+ l U 2n0 [-|S_-_ -c ] 
d n=± ' 
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The preceding equations were programmed for a digital computer and 
evaluated for various ratios of p/a The results are tabulated m 

t * 

Appendix B and agree very well with the results of an approximation 
scheme previously developed by the author (Lee, 1970) The first non- 
dimensionalization of the force is included because it follows directly 
from the equations The second is independent of the size of the masses 
V m this instance denotes the volume of one hemisphere 

For the case of the rotation (p~2a) , the components of force were 
evaluated and then transformed back to the x,y,z coordinate system by 
the matrix multiplication 


F(x,y,z) = 


cos I -sin I 




L 0 sm I 


cos I 


The results are plotted m Figure 5-7 and tabulated m Appendix B 
Note that due to symmetry the x component is always zero and 


F y (-I) = -F y (l), F a (-I) = F 2 (i), and 


Vi + 1} = Vi - I} 
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Inclination, in Degrees 


Fiqure 5-4 Normalized Components of Force as a Function of Angle of Rotation 
whenp=2a?_ 




CHAPTER 6 


CONCLUSIONS AND RECOMMENDATIONS 

As is demonstrated in chapter 5 with the example of two hemispheres, 
the force between two bodies can be efficiently and accurately calcu- 
lated after performing the necessary integrations analytically 
Except for relatively simple cases, however, the integrations still 
involve a considerable amount of work The method of this thesis 
therefore appears to be best suited for a particular problem, where a 
high degree of accuracy is required, as opposed to applying it to 
situations with widely varying geometrical parameters 

It should be noted that only potential functions expressed in 
spherical harmonics have been considered here, but analogous procedures 
may be applicable to functions expressed m spheroidal or ellipsoidal 
harmonics This is a question for further investigation 
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APPENDIX A 

GENERATION OF THE SOLID SPHERICAL 
HARMONICS <;L) H^{r,cos 6,<j>) AND 
(2) tf"(r,Cos 6,«f) 


Because all solid spherical harmonics satisfy Laplace * s equation 


= 0 , 


(A-l) 


it is useful to examine expressions which satisfy equation (A-l) m 
order to generate a particular set of solid spherical harmonics In 
this case, solutions of the form 


= (ax+by+cz) 


* - 0 , 1 , 2 , 


are examined 


Applying (A-l) , it is apparent that is a solution whenever the 
constants a,b,c satisfy 


2,2, 2 n 

a +b +c =0 


(A-2) 


By substitution into (A-2) , it can be easily verified that the following 
choice of a,b,c satisfies (A-2) 


a = l-t‘ 


b = -i(l+t z ) 


c = -2t 
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where t is a free parameter Hence 


= [ (x~iy)~ (x+iy) t 2 -2zt] * 


(A-3) 


Consistent with Figure (2-2), x,y,z are related to r,0,<j> by 


x = r sm 0 cos $ 


y = r sm 0 sm <J) 


z = r cos 0 


With these substitutions, equation (A-2) becomes 


y ^ = [r sm 0e 1 ^ ) -r sm 0e 1< ^t^-2rt cos 8]^ 


(A-4) 


With the definition of three new variables (Courant, Hilbert, 1966, 


Appendix to Ch 7) 


s = te 

s _ e l(03+7T/2) 
\1 = COS 0 


t ^ e l(w+TT/2-({>) 


(A-5) 


equation (A-4) is 


T* - (jj— ) ^ [i (p 2 -l) 1/2 (l+e 21111 ) -2xpe la> ] ^ , 


or 


= (~2rt) f & (o>) , 


(A-6) 
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where 


f A (w) = [V - cos U)(p 2 -I) 1/2 ] J1 . 


(A- 7) 


Since f^fw) is an even function of w, and for £ integer, (A-7) 
can be expanded in the following Fourier series 

a Q £ 

f *<“> - ~ + I a ra cos mw ' 

m=l 

where 

l c ¥ 

a " — / f rt (w) dto 

-IT 

i r 77 

a = — / f 0 (w)cos mtodw m = 1,2,3, , £ 

m ^ •'-Tf 


By Laplace's first and second integrals for li^> 0, n = 0,1,2, 
(MacRobert, 1947, pp. 125-129), 


r im- 


“id * i * 

i irn'P (u) 

( ji^-1) J_/ "cos w] n cos mw da) = 1 




(ntm) 1 


(A-8) 


and 


/' t »- 
J 0 


in i • t \ 

2 1/2 - n -l 1 ^{n-m'P (y) 

(y -1) x/ cos wj cos mo) do) = — t (A-9) 


These are the two integral equations which are necessary for the 
Fourier expansion of f^tw) when £ = n in the first case, and £ = -n-l 
in the second case 

The result of applying equation (A- 8) to (A-7) (£= n) is 

(1)f n (u) = V v)+2 I 1_m u+ST t p ?>> cos ™> ( A “ 10 ) 
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The range of subscripts can be extended to include negative values 
(-n ^ m ^ n) , noting that (Jackson, 1962, p 65) 


P“ m (uY = (-3-) m (n-in) 1 p m ( . 
n (n+m) ' p n lM) ' 


and thus ^f n (w) is an even function with respect to m. 
since 


sm mu = -si n(-xnw) / 


and 


p>) = p n (P), 


equation (A- 6) can be written as 


a), 


n 


(_2rt)I1 Jn P " (COS 0)eimU) 


But 


te 1 ^ = e^e 1 */ 2 


^linw = 


So, finally, 


t(x-iy) - (x+iy)t 2 -2zt) n 


" l " J ^\,n* r,ws 


Further , 


(A-ll) 
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where 


(1 ’ H n,m (r ' oos »•*> = ( ~ 21 "n;m!^ "’e 1 ”* »- 12 > 


For the set of solid spherical harmonics involving powers of r , 
l - -n-1 , 

and equation (A- 9) yields the coefficients m the Fourier expansions, 
such that 


l]x- (y^-1) ^^cos to] 


-n-1 


n m r __ \ i 

= P (y)+2 f - (n 7 m) P m (y ) cos 
n H S n' n 

m-1 


mto 


With treatment analogous to that for the first case, it is found that 


[ (x-iy) - (x+iy) t^-2zt] n ^ 


= t ~n-l V (2) 


l (r ,COS 6 r 4> ) t m , 

m = -n ' 


(A-13) 


where 


(2) H (r,cos 0 , <j>) = (n r n 1 P™ (cos 0)e lin ^ (A-14) 

11 / 111 /_2) n " t *" L n r n 
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APPENDIX B 

NUMERICAL RESULTS OF THE CALCULATION 
OF THE FORCE BETWEEN TWO HEMISPHERES 

B 1 Separation Along the Polar Axis 

Normalized Normalized 

Separation Force 


p/a 

F z /GMp d a 

V Gp a vV3 

0 007937 

-1 554 

-1 214 

0 01 

-1 549 

-1.211 

0 01587 

-I 537 

-1 201 

0 02381 

-1 521 

-1 189 

0 05002 

-1 468 

-1 148 

0 1000 

-1 376 

-1 075 

0 2 

-1 214 

-0 9485 

0 3996 

-0 963 

-0 753 

0 5 

-0 87 

-0 68 

1 998 

-0 2473 

-0 1933 

10 0 

-0 01799 

-0 01406 
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Angle of Normalized Components 

Rotation of Force 




F 

F z 

F 

y 

p z 

I (xn deg ) 

y . 

GMP d a 

z 

GMp d a 



0 

0 

0 0 

-0 2470 

0 0 

-0 1930 

22 

5 

0 01698 

-0 2529 

0.01327 

-0 1977 

45 

0 

0 03314 

-0 2711 

0 02590 

-0 2119 

67 

5 

0 04597 

-0 3019 

0.03593 

-0 2359 

90 

0 

0 05110 

-0 3416 

0 03994 

-0.2670 

112 

5 

0 04597 

-0 3814 

0 03593 

-0.2981 

135 

0 

0 03314 

-0 4121 

0 02590 

-0 3221 

157 

5 

0* 01698 

-0 4304 

0 01327 

-0 3364 

180 

0 

0 0 

-0 4363 

0 0 

-0 3410 
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NEW TECHNOLOGY APPENDIX 


This report is published in the belief that it constitutes 
an improvement in the state of the art In particular/ pages 
9 through 82 are referenced 



